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Direct simulations (ILES) of turbulent, separated flow over the wall-mounted hump 
configuration is conducted to investigate the physics of separated flows. A chord-based 
Reynolds number of Re c = 47,500 is set up, with a turbulent inflow of Reg = 1,400 (6/c 
= 3%). FDL3DI, a code that solves the compressible Navier-Stokes equations using high- 
order compact-difference scheme and filter, with the standard recycling/rescaling method of 
turbulence generation, is used. Two different configurations of the upper-wall are analyzed, 
and results are compared with both a higher Re c (= 936,000, Reg = 7,200, 6/c = 0.77%) 
experiment for major flow features, and RANS (k-cu SST) results. A lower Re c allows 
for DNS-like mesh resolution, and an adequately wide span. Both ILES and RANS show 
delayed reattachment compared to experiment, and significantly higher skin friction in the 
forebody of the hump, as expected. The upper-wall shape influences the C v distribution 
only. Results from this study are being used to setup higher Re c (lower 6/c) ILES. 


I. Introduction 

Accurate prediction of turbulent flows at high Reynolds numbers continues to be a significant problem in 
fluid dynamics. When turbulent boundary layers separate, predictions become all the more challenging since 
our understanding of the underlying physics is still lacking. Currently, most turbulent flows in engineering 
applications are simulated using various models that represent the understood physics. There is still no 
universal model that can be used to compute every turbulent flow scenario. 

Various workshops over the past decade have brought together CFD practitioners to assess current ca- 
pabilities of different types of turbulent flow solution methodologies, and use benchmarked experiments to 
identify potential areas of improvement. In 2004, the Langley Research Center Workshop on CFD Validation 
of Synthetic Jets and Turbulent Separation Control 1,2,3 (CFDVAL2004) considered three flows of interest 
that had been studied experimentally, and had extensive statistical databases of mean and turbulent quan- 
tities. Case 3 at the workshop was flow over the wall-mounted hump configuration, which consisted of a 
Glauert-Goldschmied-type airfoil mounted between two glass end-plates. It is a nominally two-dimensional 
experiment, although there are three-dimensional side-wall effects near the end-plates. This is based on ear- 
lier experiments of Seifert and Pack, who had studied flows over a similar configuration at higher Reynolds 
numbers (Re c = 2.4-2.6xl0 6 , based on chord-length). The baseline experiment that was used as the refer- 
ence in the workshop did not involve flow control, and was conducted by Greenblatt et al. * 1 It had an Re c 
of 936,000, with a turbulent inflow Reg = 7,200 at = 0.1. Both leading and trailing edges were faired 
smoothly with a wind tunnel splitter plate, whereas the experiments of Seifert and Pack installed the airfoil 
on the wind tunnel floor. This baseline experiment without flow-control showed that flow separates near 65% 
of the chord, and reattaches downstream past the hump. Other than C p and Cf distributions, 2D/stereo 
PIV data of velocity profiles and turbulence quantities along the tunnel center plane were documented. 

At the workshop, this case was analyzed extensively. Numerous RANS (Reynolds-averaged Navier- 
Stokes) studies were conducted and all failed to reproduce the many characteristics of this flow. Iaccarino et 
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al. used CIRA a -UZEN b and the commercial code Fluent to test the Spallart-Allmaras and k-w turbulence 
models, and could not predict this baseline flow with satisfactory agreement. Capizzano et al. used a 
modified version of CIRA-UZEN, and tested the Menter-SST k-ui model. They found that the Reynolds 
stresses downstream of the point of separation were under-predicted, and reattachment therefore delayed. 
Cui and Agarwal used an unsteady RANS solver (WIND) and the SST model, and found the combination 
to severely over-predict the surface pressure. This also resulted in a longer separation bubble. Each of these 
results used a second-order structured-grid methodology to compute the flow. Balakumar used a fifth-order 
WENO (weighted essentially non-oscillatory) , with third-order TVD (total variation diminishing) Runge- 
Kutta, and the SST model. Morgan et al. used second- and fourth-order compact-difference scheme and 
the k-e turbulence model for the case. Both resulted in a similar under-prediction of the Reynolds stresses, 
and a longer bubble. Bettinni and Cravero, 1 He et al. and Madugundi et al. used Fluent with different 
turbulence models. Bettinni and Cravero used the RNG k-e model, and found it to significantly under-predict 
surface pressure over the attached portion of the hump, and a very rapid rise in pressure in the separated 
region. He et al. used the k-e, S-A, and SST models, and made significantly better pressure predictions using 
the same software. S-A and SST models still over-predicted the separation bubble size, and the k-e model 
predicted later separation. Madugundi et al. surprisingly under-predicted the size of the separation bubble 
using the same software, with a number of turbulence models. The all-pervading observation with RANS 
models is the inconsistency of the size of the separation bubble, attributed to eddy-viscosity and Reynolds 
stresses being under-predicted within the separation region. Rumsey 1 ! used CFL3D, a URANS code, with 
three different turbulence models, and showed that doubling eddy-viscosity in a region near the separation 
bubble did indeed result in earlier reattachment, more comparable to experiment. 

Some LES and hybrid RANS/LES methods made much better predictions of the flow features. Biswas' 
used an incompressible fifth-order WENO code with upwind-biased finite differencing to solve the LES 
equations. In this, sub-grid scale turbulent stresses were closed using a dynamic turbulent kinetic energy 
transport model. A spanwise width of 0.2c (c - chord length) was set with periodic boundaries. Very good 
agreement of surface pressure and skin friction were seen with experiment. Krishnan et al. conducted 
RANS simulations (using S-A and SST), as well as DES of the hump flow problem, and showed that DES 
could improve the bubble size prediction, reducing it by 20% compared to RANS. The velocity profiles 
compared well with experiment, but the Reynolds shear stress <uv> was slightly over-predicted. Their DES 
was setup with a periodic span of 0.121c. Saric et al. solved the incompressible RANS, LES and DES 
equations using a cell-centered collocated finite-volume scheme, with the SIMPLE algorithm. They used 
the S-A model for RANS and for the unresolved scales in DES (Az + =150), and the Smagorinsky sub-grid 
scale model for LES with a grid spacing of Aa; + = 80, A y + = 1, A z + = 50 at the wall. The spanwise widths 
were 0.2c for DES and 0.152c for LES (which was found to be insufficient to ensure spanwise decorrelation), 
with periodic boundary conditions. In spite of the coarseness, LES showed very good comparisons with 
experiment in terms of shear stresses and bubble length. DES performed similarly well, and the S-A RANS 
simulation once again over-predicted the bubble length. You et al. performed incompressible LES with an 
energy- conservative second-order central difference scheme on a staggered mesh, where aliasing errors were 
controlled by enforcing kinetic energy conservation. This was done with the dynamic Smagorinsky sub-grid 
scale model. A recycling inflow procedure was used on a mesh with a spanwise extent of 0.2c, and a wall 
resolution of Ax + < 50, Ay + < 0.8, Az + < 25. Surface pressure, mean velocity profiles and Reynolds stresses 
compared very well with experiment. Skin friction in the fore-body of the airfoil was over-predicted. But, the 
size of the separation bubble agreed well with experiment. Franck and Colonius ' conducted compressible 
LES and ILES with a non-density weighted low-pass filtered method, with a sixth-order Pade scheme in the 
wall-normal direction, fourth-order explicit treatment in the streamwise direction, and a Fourier spectral 
method in the spanwise direction. They used a filtering optimization parameter of a/ = 0.45. The mesh 
resolution was Ax + < 94, Ay + < 8.7, A z + < 31, with a spanwise extent of 0.2c. A much finer resolution was 
imposed over the hump, and within the separation region. Overall, the velocities, Reynolds stresses and 
the bubble length agreed well with experiment. They compared results using both constant coefficients and 
dynamic Smagorinsky, as well as using no sub-grid scale model, and found very little difference. They also 
analyzed results at higher Mach numbers: M^ = 0.1 and 0.3 were identical, whereas M^ = 0.6 showed 
higher acceleration over the hump, and a longer separation bubble. 

There were other hybrid RANS/LES that did not improve on the RANS results during the workshop. 
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Israel et al. used the blended RANS/LES method flow simulation methodology (FSM) in a second- 
order code based on flux difference splitting and dual time-stepping. The underlying RANS model used 
was an explicit algebraic stress model (EASM) based on k -e. With a span of 0.17c, the time-averaged 
results showed a reattachment point too far downstream, similar to RANS. Hiller and Seitz 21 used the 
scale adaptive simulation (SAS) blended RANS/LES model on a domain with 0.56c span. Their limited 
quantitative comparisons showed pressure distribution similar to the k -ui RANS model, with much higher 
pressure within the separation region than experiment. 

Morgan et al. conducted compressible ILES with a fourth-order compact-difference scheme, and used 
a high-order non-dispersive Pade-type filter instead of a sub-grid scale model. Their method was similar 
to that of Franck and Colonius, with ctf = 0.4 in a sixth-order filter. To attain a reasonable resolution 
(A.t + = 60,Ay + = 1,A z + = 20), they scaled down Re c to 200,000 on a mesh with a spanwise width of 
1.39c. This corresponded to a momentum thickness to chord ratio ( 9/c ) of 3.6%. The mean velocity profiles 
compared well except near the point of reattachment, and the suction pressure was higher than experiment, 
as was skin friction in the fore-body of the hump. These were attributed to the lower Reynolds number used 
in the simulation. Upstream of the separation point, Reynolds stresses were significantly over-predicted, and 
it should be noted that for this lower Reynolds number, the separation point was further upstream compared 
to experiment. 

Postl and Fasel conducted the only known DNS for this case (using the incompressible vorticity- velocity 
form of the Navier-Stokes equations), with a mesh resolution that resulted in the simulation being under- 
resolved (Aa; + = 27 — 92, A y + = 1.2, Az + = 17), and a span of 0.142c. In addition to being under-resolved, 
a thinner laminar inflow profile with Reg = 4,000 at x/c = -0.5 was used. This is equivalent to a Reg = 
2,400 at x/c = -2.14 ( 9/c = 0.26%). This was tripped near the inflow using a time-harmonic forcing term. 
Even with that setup, the separation point and mean velocity profile before separation matched well with 
experiment. But the size of the separation bubble and reattachment point were both over-predicted by about 
20%. This was attributed to the relative coarseness of the mesh, which failed to resolve the smallest scales 
adequately. Further, lack of information about freestream turbulence intensity added to these discrepancies, 
particularly in the range of frequencies where the separated shear layer is hydrodynamically unstable. 

Besides CFDVAL2004, more recent workshops have also included this test case to further assess the 
state-of-the-art in CFD methods. Both the 11th and the 12th ERCOFTAC c /IAHR d Workshops on Refined 
Turbulence Modelling had results of flow over the hump-flow configuration. ' ’ This test case has since 
been added to the ERCOFTAC Database (Classic Collection) as Case C.83®. It was also included in the 
ECCOMAS CFD 2010 Conference as a part of the V&V I: Code Verification and Solution Verification & 
Validation Workshop f , and the Advanced Turbulence Simulation for Aerodynamic Application Challenges 
(ATAAC) project 8 . Further interest in this test case has even led to lattice Boltzmann simulations of the 
configuration, coupled with a very large eddy simulation (VLES) turbulence model. Noelting et al. used 
a 19-state lattice Boltzmann model (D3Q19), with the collision integral in the Boltzmann equation being 
approximated using the Blratnagar-Gross-Krook (BGK) form. Turbulence modeling was directly included 
in the lattice Boltzmann equations by modifying the relaxation time, to give an extended relaxation time 
scale. This extended time was then calculated using a variant of the RNG k-e model. This resulted in 
the emergence of significantly more complex physics from the RNG k-e equations, than when using the 
Navier-Stokes equations. Results were similar to DES and LES in terms of surface pressure and bubble 
length. 

Other LES and hybrid RANS/LES methods have also been studied in more detail since the workshop, and 
many efforts have shown very good progress and excellent comparisons with experiment, ’ ’ ’ especially 
with DDES. ' But, the significant shortcomings of RANS models are yet to be addressed. Yoshio and 
Abe ’ and Cappelli and Mansour recently documented these failures, including the under-prediction of the 
Reynolds stresses within the separation region, a farther reattachment location, and a slower recovery of the 
boundary layer after reattachment. 

Seifert and Pack 1 noted that for the given fully turbulent inflow conditions, the characteristics of separa- 
tion and the bubble length were largely insensitive to the inflow Reg and Re c . This was especially true since 
laminar to turbulent transition was eliminated in the domain. But the ratio of the thickness of the incoming 
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boundary layer to the hump height (S/h) played a significant part. In other words, matching experimental 
9/c is critical in replicating the conditions for similar flow separation. Post and Fasel, on the other hand, 
attributed differences between their DNS and experiment to the coarseness of their mesh, and a narrow 
span. Given the computational costs of setting up DNS at the experimental conditions (estimated to be 
upwards of 6 billion grid points with the current setup), a preliminary test case of lower Reg:Re c flow over 
the wall-mounted hump is considered. In compromising on the 9/c ratio, a mesh resolution comparable to 
DNS is realizable, as is a wide-enough span. This is conducted as a preliminary step towards setting up 
subsequent simulations at experimental 9/c , in addition to lending an extra data-point, considering the lack 
of high quality DNS for this configuration. RANS simulations using the k-w SST model at this Reg:Re c 
combination are also conducted. Comparisons of pressure and skin friction, as well as velocity and Reynolds 
stress profiles are made with the available experimental data with a lower 9/c. 

This paper is organized as follows: Section II briefly explains the numerical methodologies used for 
the ILES. Section III presents the results, which include the turbulent inflow generation using the recy- 
cling/rescaling method, as well as flow over the wall-mounted hump geometry. Qualitative and quantitative 
comparisons are made with experiment, as well as with RANS. Conclusions are drawn in Sec. IV. 

II. Numerical Method 

II. A. Governing Equations, Compact Difference Scheme and Filter 

The finite-difference code (a variant of AFRL’s FDL3DI l ’’ i ' ) used in this study solves the compressible 
Navier-Stokes equations for an ideal gas. The sixth-order compact-difference scheme of Lele 8 is used to 
solve the governing equations in transformed curvilinear coordinates. ’ Fourth-order explicit Runge-Kutta 
(RK-4) is used for time integration. During each time step, an eighth-order low-pass spatial filtering scheme 
is applied to the conservative variables to ensure stability, along with second- and fourth-order near-boundary 
formulations. A filtering optimization parameter of a/ = 0.495 is set. This results in lesser dissipation 
compared to the setups used by Morgan et al. ; («/ = 0.4) and Franck and Colonius 1 (a/ = 0.45). To 
ensure stability and accuracy of time integration, a Courant-Friedrichs-Lewy (CFL) number in the vicinity 
of 0.4, based on the acoustic speeds, is used in the RK-4 scheme. 

This code was previously designed to solve for supersonic flows, specifically those involving shock-wave 
and boundary layer interaction problems. Therefore, the localized artificial diffusivity (LAD) method of 
shock-capturing was previously installed by Kawai et al., 42, and was used in a detailed study of the flow 
physics of such interactions by Morgan et al. This method adds artificial coefficients to the fluid transport 
terms, and takes effect in the vicinity of shocks. Since the flow regime under consideration is subsonic, LAD 
is not included in the current study. Thus, other than the low-pass filter, no additional sub-grid scale model 
is used. In that sense, this effort can be considered implicit-LES, or more precisely “under-resolved DNS.” 

II. B. Recycling/Rescaling Turbulence Inflow Method 

For generating the turbulent inflow, the standard recycling/rescaling method of Lund et al. is used in 
the code. Urbin and Knight extended this method to compressible flows for solving a Mach 3 turbulent 
boundary layer. This was implemented as a part of the effort to study the interaction of an oblique shock 
wave and a turbulent boundary layer by Morgan et al. The approach extracts the instantaneous velocity 
profile from a plane downstream of the inflow, rescales and reintroduces it as the inflow. The schematic 
of the procedure is shown in Fig. 1. It has been shown to produce a realistic turbulent boundary layer, 
yielding statistics that are in good agreement with experiment and DNS. To further prevent contamination 
of the solution with spurious spatio-temporal correlations generated by this reintroduction procedure, a non- 
constant reflection of the recycled turbulence plane is applied at randomly distributed time intervals. This 
method, called recycling/rescaling with dynamic reflection (RR+DR), has been used to solve subsonic flows 
here. 
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III. Results 


III. A. Computational Setup 

Details of the wall-mounted hump simulation are explained here. This is followed by validation of the inflow 
profile, and the comparison of results of the ILES with RANS and experiment. 

As mentioned before, the presented ILES corresponds to a higher 9/c on a mesh with DNS-like resolution. 
For the inflow to the hump domain, a turbulent velocity profile at = 0.3 and Reg = 1,400 was chosen. 
This is in contrast to the experimental inflow of = 0.1 and Reg = 7, 200. Franck and Colonius showed 
that the flow over the hump is insensitive to differences in freestream Mach number, from 0.1 to 0.3. Reg was 
so chosen to allow for the relatively inexpensive validation of the inflow profile with the spectral DNS data 
of Spalart. 1 With these conditions, a domain of length 156r, height 35r and span 15r was used, where 6 r 
is the desired boundary layer thickness at the velocity capture-plane. The height of the domain was chosen 
based on the validation study of Lund et al. for the same Re#. The spanwise width was identified in 
a parametric study that investigated the minimum width required for turbulent boundary layers to match 
reference DNS. 



Figure 2. Perspective view of the schematic of the computational setup 


The dimensions of the hump domain were defined based on the CFDVAL2004 workshop. The inflow 
plane was located 2.14 chord-lengths (c) upstream of the leading edge of the hump, and the exit plane 4.0c 
downstream, respectively. The height of the domain from the splitter-plate to the upper wall was 0.90905c. 
The chord-based Reynolds number was computed based on the constraints set by the dimensions of the inflow 
and hump domains. Thus, 3<5r = 0.90905c (Fig. 2). The theoretical <5r for this inflow, based on Prandtl’s 
one-seventh power law, was computed to be 0.002061 m, resulting in a chord of 0.0068 m. The experimental 
chord length was 0.42 m. With a freestream Mach number of 0.3, the ILES was setup for an Re c ~ 47, 500, 
which is one-twentieth of the experimental Re c (= 936,000). Thus, the 9/c ratio was computed to be 3%, 
compared to an experimental 9/c of 0.77%. It should be mentioned here that the coarse DNS of Postl and 
Fasel had a 9/c of 0.26%, whereas the lower Re c ILES of Morgan et al. had their 9/c = 3.6%. With such 
significant differences in flow conditions, major flow features and trends in pressure and skin friction were 
the main points of comparison in this study, as were the mean velocity profiles and Reynolds stresses. 
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III.B. Turbulent Inflow 


The turbulent inflow profile with R eg = 1, 400 was generated simultaneously as an auxiliary simulation using 
the recycling/rescaling method. ’ ’ The inflow capture-plane was located about 12<5p{ downstream of the 

inlet, and was also the location of the target Reg (Fig. 1). This mesh consisted of 500 x 100 x 105 points, 
and was uniformly spaced along the streamwise and spanwise directions. A hyperbolic tangent stretching 
was set in the wall normal direction. With these, the mesh resolution at the wall was comparable to DNS, 
with Ax + = 19.4, A?/ + =0.65, A;z + = 6. This was validated by comparing the mean velocity profile here with 
the spectral DNS data of Spalart (Fig. 3). 



Figure 3. Mean velocity profile in wall units, at x/c = -2.14, until = 1000, compared with DNS results for fiat plate 
of Spalart 


There is very good agreement through the viscous sublayer, and the log layer. The slightly higher 
freestream velocity at the edge of the boundary layer is attributed to the upper-wall boundary condition. 
Whereas the spectral DNS corresponds to a zero pressure-gradient flat plate boundary layer, the upper-wall 
in the inflow domain is set to a symmetry boundary. This results in a slight acceleration of the freestream. 
The Reynolds stresses at the same 125 location are compared with reference spectral DNS in Fig. 4. 

Figure 4 (a) compares the Reynolds stresses throughout the boundary layer. Differences are evident 
within the buffer zone, whilst the rest of the inflow profile compares well with spectral DNS. The over- 
prediction, also evident in the mean velocity profile, is discussed in a separate paper by the authors, 1 and 
is once again attributed to differences in the upper-wall boundary condition. 

III.C. Wall-mounted Hump 

With the validated inflow profile, ILES of the hump flow was carried out on two separate meshes, shown 
in Fig. 5. These meshes differed in the shape of the upper-wall. During the CFDLVAL2004 workshop, the 
top-wall shape was so modified to account for side-wall blockage effects, which improved the agreement of 
the pressure distribution over the hump with experiment. At this higher 9/c , both upper- wall profiles were 
tested to document differences at the current flow conditions, “structl” corresponds to the flat upper-wall 
configuration, and struct4 has the modified shape. Each mesh consisted of 827 x 100 x 105 points, with a 
similar DNS-like resolution at the wall (Aa; + < 19.4, Ay + < 0.65, Az + < 6). It should be noted here that the 
coarse DNS of the higher Re c flow by Postl and Fasel used a mesh with Arr + = 27 — 92, A y + = 1.2, A z + = 17, 
on a domain of span 0.142c. The current setup had a span of 0.3c. This setup, therefore, address both 
drawbacks of the earlier coarse DNS, despite being for a higher 9/c. 

For the remaining boundary conditions, a freestream pressure outflow was set at the outlet. No-slip, 
adiabatic walls were used for the bottom and top walls. Whereas most of the simulations during the workshop 
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(a) Across the boundary layer 


(b) In the viscous near-wall region 


Figure 4. Reynolds stresses normalized by friction velocity, at x/c — -2.14, compared with DNS results of Spalart 


used a symmetry or slip-wall boundary at the upper-wall, this was so chosen to be a better representation 
of the experimental setup. Along the span, periodic boundary conditions were set. 

The non-dimensionlized mean velocity contours are shown in Fig. 6. Figure 6 (a) shows the experimental 
2D PIV results for the Re c = 936,000 case, and Figs. 6 (b) and (c) represent the ILES results for Re c = 
47,500 using meshes “struct 1” and “struct4,” respectively. The first observation that stands out is the longer 
separation bubble at the higher 0/c (3%). There is little difference between the two meshes, and both predict 
the reattachment point downstream of that for the experiment with 9/c = 0.77%. In addition, the point of 
separation is further upstream. Both these observations conform with those seen in the ILES of Morgan et 
al., whose 0/c was a comparable 3.6%. The recirculation bubble is also evident in Fig. 7, which represents 
the streamlines of the flow, shown for the experiment (a), and ILES ((b) and (c)). The earlier separation 
and delayed reattachment are both clearly apparent here. 

For a quantitive comparison of the flow features, pressure and skin friction distributions along the wall 
were compared with experiment, and RANS results. OpenFOAM h (Open source Field Operation And 
Manipulation), used for the RANS simulations, is a free, open source software suite written in C++ that 
allows for the development of customized numerical solvers and pre/post-processing tools for solving fluid 
dynamics or other continuum mechanics problems. Solving turbulent flows requires a Navier-Stokes equations 
solver that includes a turbulence model to close the Reynolds stress terms. For these RANS simulations, 
simpleFoam , a steady-state solver for incompressible turbulent flows that employs the Semi-Implicit Method 
for Pressure Linked Equations (SIMPLE) algorithm was used with the k — uj SST model. Suitably scaled 
down versions of CFDVAL2004 workshop meshes (scale = 1/20) were used, with a steady freestream inflow 
set at x/c = -3 (instead of x/c = -6.39 used at the workshop) to account for the reduced Re# at the inlet. 
Simulations with both upper-wall shapes were tested, with both slip and no-slip upper-wall boundaries. 

Comparisons of C p and Cf are shown in Fig. 8. In Fig. 8 (a), pressure over the fore-body of the hump are 
significantly higher for ILES and RANS than the Re c = 936,000 experiment, and lower within the separation 
bubble. The higher pressure over the fore-body was also seen in the ILES of Morgan et ah, but C p within 
the separation bubble is counter to their observation. Whereas their ILES showed higher pressure within the 
separation bubble compared to experiment, the current ILES showed an even lower C p within the bubble, 
as did RANS. Also, in agreement to the cases in the CFDVAL2004 workshop, the modified upper-wall did 
result in a lower C p over the fore-body section of the hump. 

“www.openfoam.com 
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(a) struct 1: default flat upper- wall 



(b) struct4: modified upper-wall, contoured to match experimental pressure profile 

Figure 5. Meshes used for ILES: 827x100x105 points 


Figure 8 (b) shows the variation of skin friction for the experiment, ILES and RANS. Once again, similar 
to the observation of Morgan et ah, the higher 9/c results in a significantly higher Cf in the fore-body of the 
hump. Whereas their peak Cf was 0.01, ILES of the current flow conditions had a peak Cf of 0.015. Skin 
friction with RANS peaked at 0.018. Based on this plot, the point of reattachnrent for ILES was estimated 
to be at x/c = 1.25, and ss 1.4 for RANS. This is in contrast to the experimental result of x/c ss 1.1 for 9 /c 
= 0.77%. There is, also, virtually no difference between both ILES meshes with regards to skin friction. 

Mean streamwise and wall normal components of the flow velocity at different locations downstream of 
separation are compared in Fig. 9. Experimental data at x/c = 1.0, 1.1, 1.2 and 1.3 are contrasted with 
ILES for both upper-wall definitions. Marginal differences between the upper-wall shapes are noticed only 
in the freestream. Figure 9 (a) shows the velocity profiles at x/c = 1.0. ILES shows a stronger negative 
streamwise component compared to experiment. This trend continues at other locations downstream. Based 
on the experiment, for 9/c = 0.77%, the flow reattaches at x/c ss 1.1, which is evident in Fig. 9 (b). For 9/c 
= 3%, the flow continues to be separated until x/c ~ 1.3, where recovery starts, as in Fig. 9 (d). The point 
of reattachnrent based on these plots should lie between x/c = 1.2 and 1.3, which agrees with the value x/c 
= 1.25 from the Cf plot in Fig. 8 (b). 

In Fig. 10, mean Reynolds stress components <uu > , <vv> and <uv> are compared at the same x/c 
locations with experiment. <uu> and <vv> compare well at x/c = 1.0 until y/c = 0.6, above which the 
thicker inflow profile ( 9/c = 3%) tends to exhibit higher magnitudes than the thinner profile ( 9/c = 0.77%). 
A similar trend is observed for <uu> and <uv> at the downstream locations as well, whereas the wall 
normal component <vv> is significantly higher for the thicker profile throughout the boundary layer (except 
very close to the wall) . The higher magnitude of Reynolds stresses conforms with the observation of a longer 
separation bubble. 
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(a) PIV for Re c = 936, 000 
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(b) ILES (struct 1) for Re c = 47, 500 



J i i i I i i i I i i i I 


0.8 1 1.2 1 . 

(c) ILES (struct4) Re c = 47, 500 

Figure 6. Contours of U/Um 


Considering the differences in inflow profiles between the simulations and experiment, to make a better 
estimate of recovery, the axial location of reattachment for the experiment was shifted so as to coincide with 
ILES, and the mean profiles and Reynolds stresses downstream of that point were compared. These plots 
are shown in Figs. 11 and 12. 

The overall trend of recovery shows very good agreement for both 9/c inflow profiles. The experiment 
with the thinner profile is slightly offset to the right, which is attributed to the differences in Reynolds 
numbers. But the shapes of both U and V match very well. Similarly, there is significantly improvement 
agreement of the Reynolds stresses, evident in Fig. 12. Here again, differences are attributed to the inflow 
conditions, but the trend of the stresses from ILES while progressing downstream conform with experiment 
very well. 


IV. Conclusions 

As a part of the effort to improve turbulence models for separated flows, and considering the lack of high 
quality DNS for the wall-mounted hump configuration, direct simulations of flow with thicker inflow profiles 
have been conducted. By compromising on the 9/c ratio (3% for simulations vs 0.77% for experiment), 
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(a) PIV for Re c = 936, 000 



(b) ILES: structl (Re c = 47, 500) 



(c) ILES: struct4 (Re c = 47, 500) 

Figure 7. Streamlines 


DNS-like mesh resolutions, as well as spans of adequate widths, could be afforded. 

Qualitative comparisons with experiment show an earlier point of separation for the higher 9/c simu- 
lations, and a later point of reattachment. The overall trends of pressure and skin friction variations over 
the surface agree, but the magnitudes are vastly different. The thicker inflow profiles result in significantly 
higher pressures in the fore-body of the hump, and lower pressures within the separation bubble. 

Comparisons of the mean velocity profiles along different streamwise locations downstream of reattach- 
ment once again show differences in the bubble size, as well as recovery. Whereas the experimental reattach- 
ment point ( x/c ) is 1.1 chord-lengths downstream of the leading edge of the hump for the thinner inflow, 
the simulation estimates it at x/c = 1.25 for the thicker inflow. Reynolds stresses close to the wall tend to 
agree for both inflow conditions. But considering the differences in flow profiles, direct comparisons between 
experiment and simulations are inconsequential. 

Instead, by shifting the streamwise location of reattachment to match with simulations, the trends of 
recovery, both in terms of mean velocity profiles and the Reynolds stresses, could be compared. Even though 
the magnitudes differ marginally because of the flow Reynolds numbers, the overall recovery trends in terms 
of mean velocities agree well for both inflow profiles. A similar agreement is also seen for the Reynolds 
stresses. 

Two meshes with different upper-wall shapes were used, and comparisons of mean velocity profiles and 
Reynolds stresses show near identical values downstream of the point of separation. The skin friction 
variation along the surface is also the same. There are differences, though, in the pressure distribution over 
the surface of the hump, as well as within the separation bubble. But, this clearly has very little impact on 
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Figure 8. Pressure and skin friction variation 




Figure 9. Mean velocity profiles (U/Uqo and V/Uqo) at different streamwise locations 


the characteristics of separation, be it the point of separation or reattachment, as well as recovery. 

Based on the inferences from this initial test, additional simulations with inflows comparable to exper- 
iment are being conducted. Three-dimensionality effects of the side-walls, which the modified upper-wall 
shape was intended to account for, will also be investigated. 
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